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Abstract 

Let G(n,c/n) and G r (n) be an n-node sparse random graph and a sparse random r- 
regular graph, respectively, and let T(n, r) andZ(n, c) be the sizes of the largest independent 
set in G(n, c/n) and G r (n). The asymptotic value of T(n, c)/n as n — > do, can be computed 
using the Karp-Sipser algorithm when c < e. For random cubic graphs, r = 3, it is only 
known that .432 < lim inf n I(n, 3)/n < limsup n T(n, 3) < .4591 with high probability 
(w.h.p.) as n — > oo, as shown in |FS94| and |Bol81j . respectively. 

In this paper we assume in addition that the nodes of the graph are equipped with non- 
negative weights, independently generated according to some common distribution, and 
we consider instead the maximum weight of an independent set. Surprisingly, we discover 
that for certain weight distributions, the limit lim n X(n, c)/n can be computed exactly even 
when c > e, and lim n I(n,r)/n can be computed exactly for some r > 2. For example, 
when the weights are exponentially distributed with parameter 1, lim n T(n, 2e)/n ~ .5517, 
and lim n T(n, 3)/n .6077. Our results are established using the recently developed local 
weak convergence method further reduced to a certain local optimality property exhibited 
by the models we consider. Using the developed technique we show in addition that in the 
unweighted case liminf n X(n,4)/n > .3533, which is a new lower bound. We also prove 
that in any (non-random) graph with degree 3 and large girth, the size of the maximum 
independent set is at least .3923n — o(n), improving the previous bound (7/18)n — o(n) 
in HS82 . Finally, we extend our results to maximum weight matchings in G(n,c/ri) and 
G r (n). For the case of exponential distributions, we compute the corresponding limits for 
every c > and every r > 2. 
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1 Introduction 

Two models of random graphs considered in this paper are a sparse random graph G(n,c/n) 
and a sparse random regular graph G r (n). The first is a graph on n nodes {0, 1, . . . , n — 1} = 
[n], where each potential undirected edge (i,j),0 < i < j < n — 1 is present in the graph 
with probability c/n, independently for all n(n — l)/2 edges. Here c > is a fixed constant, 
independent of n. A random r-regular graph G r (n) is obtained by fixing a constant integer r > 2 
and considering a graph selected uniformly at random from the space of all r-regular graphs on n 
nodes (graphs in which every node has degree r) . A set of nodes V in a graph G is defined to be 
an independent set if no two nodes of V are connected by an edge. Let Z(n, c) and Z(n, r) denote 
the maximum cardinality of an independent set in G(n,c/n) and G r (n) respectively. Suppose 
the nodes of a graph are equipped with some non- negative weights Wi,0 < i < n — 1 which 
are generated independently according to some common distribution F w (t) = P{W / j < t},t > 0. 
Let X w (n,c),I w (n,r) denote maximum weight of an independent set in G(n,c/n) and G r (n) 
respectively. 

A matching is a set of edges A in a graph G such that every node is incident to at most 
one edge in A. Let A4(n,c) and Ai(n, r) denote the maximum cardinality of a matching in 
G(n, c/n) and G r (n), respectively. It is known that G r (n),r > 3 has a full matching w.h.p., that 
is Ji4(n,r) = n/2 ([n/2\ for odd n) w.h.p. [JLR00J. If the edges of the graph are equipped 
with some non-negative random weights, then we consider instead the maximum weight of a 
matching A4 w (n,c) and A4 w (n, r) in graphs G(n, c/n) , G r (n) , respectively. The computation of 
I w (n,c),I w (n,r),A4 w (n,c),A4 w (n,r) in the limit as n — > oo is the main subject of the present 
paper. 

The asymptotic values of X(n, c) for c < e and A4(n, c) for all c were obtained by Karp and 
Sipser using a simple greedy type algorithm in |KS81j . Yet it is an open problem to compute 
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the corresponding limit for for independent sets for the case c > e or in random regular graphs 
or even to show that the limit exists [AldJ. 

The developments in this paper show that, surprisingly, proving the existence and the com- 
putation of the limits lim n X(n, -)/n, lim n A4(-)/n is easier in the weighted case than in the un- 
weighted case, at least for certain weight distributions. In particular, we compute the limits for 
independent sets in G r (n), r = 2, 3, 4 and G(n, c/n), c < 2e, when the node weights are exponen- 
tially distributed, and we compute the limits for matchings in G r (n) and G(n,c/n) for all r, c, 
when the edge weights are exponentially distributed. It was shown by the first author |Gam02| 
that the limit lim n JH(r, c)/n exists for every weight distribution with bounded support, though 
the non-constructive methods employed prevented the computation of the limits. 

Our method of proof is based on a powerful local weak convergence method developed by 
Aldous [Ald92j, [AldOlJ, Aldous and Steele AS03J, Steele [Ste02 j; further empowered by a certain 
local optimality observation derived in this paper. Local weak convergence is a recursion technique 
based on fixed points of distributional equations, which allows one to compute limits of some 
random combinatorial structures, see Aldous and Bandyopadhyay [AB for a recent survey on 
applications of distributional equations and Aldous and Steele |AS03| for a survey on the local 
weak convergence method. In particular, the method is used to compute maximum weight 
matching on a random tree, when the weights are exponentially distributed. The tree structure 
was essential in |AS03| for certain computations and the approach does not extend directly to 
graphs like G(n,c/n) with c > 1, where the convenience of a tree structure is lost due to the 
presence of a giant component. It was conjectured in |AS03j that a some long-range independence 
property might be helpful to deal with this difficulty The present paper partially answers this 
qualitative conjecture in a positive way. We introduce a certain operator T acting on the space 
of distribution functions. We prove a certain local optimality property stating that, for example, 
for independent sets, whether a given node i belongs to the maximum weight independent set is 
asymptotically independent from the portion of the graph outside a constant size neighborhood 
of i, iff T 2 has a unique fixed point distribution. Moreover, when T 2 does have the unique fixed 
point, the size of the extremal object (say maximum weight independent set) can be derived 
from a fixed point of an operator T. The computations of fixed points is tedious, but simple in 
principle and the groundwork for that was already done in |AS03j . We hope that the long-range 
independence holds in other random combinatorial structures as well. 

The issue of long-range independence of random combinatorial objects is addressed in a some- 
what different, statistical physics context in Mossel [Mos03 , Brightwell and Winkler BW03J, 
Rozikov and Suhov |RS03j . Martin [Mar03j, where independent sets (hard-core model) are con- 
sidered on infinite regular trees, weighted by the Gibbs measure. The long-range interaction 
between nodes at a large distance is investigated with respect to this measure using the notion of 
reconstruction. It would be interesting to investigate the connections between the two models. 

In a different setting Talagrand [Tal03j proves a certain long-range independence property for 
the random assignment problem, where the usual min-weight matching is replaced by a partition 
function on the space of feasible matchings. He uses a rigorous mathematical version of the 
cavity method, which originated in physics, to prove that the spins (edges of the matching) are 
asymptotically independent as the size of the problem increases. The particular form of the long- 
range independence is similar to the one we obtain, refer to Theorem El below, and in fact the 
cavity method, which is based on "knocking" out certain spins from the system and analyzing 
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the relative change of the size of the extremal object, has some similarity with the local weak 
convergence method, which is also based on considering extremal objects (say independent sets) 
with one or several nodes excluded. It seems worth investigating whether there is a formal 
connection between the two methods. Finally, we refer the reader to Hartmann and Weigt 
HWOlJ who derive the same result as Karp and Sipser for independent sets using non-rigorous 
arguments from statistical physics. 

The rest of the paper is organized as follows. In the following section we describe some prior 
results on maximum independent sets and matchings in random graphs. Our main theorems are 
given in Section El The operator T, fixed point equations and long-range independence issues 
are discussed in Section HJ The main results are proven in Section Some conclusions are in 
Section |H1 

We finish this section with some notational conventions. Exp(/i), Pois(A), Be(p) denote re- 
spectively exponential, Poisson and Bernoulli distributions with parameters /i, A > 0, < p < 1. 



2 Prior work and open questions 

It is known and simple to prove that X(n,r) = I(n,c) = Q(n) w.h.p. for any constants r > 
l,c > 0. Moreover, it is known that, w.h.p., 61og(3/2) — 2 = .432... < liminf n X(n, 3)/n < 
limsup n X(n, 3)/n < .4591. The lower bound is due to Frieze and Suen [FSf) 1 , and the upper 



bound is due to Bollobas [B0I8IJ. The upper bound is generalized for any r > and uses a very 
ingenious construction of random regular graphs via matching and random grouping, |Bol80j . 
|.)LK(K)| . 

It is natural to expect that the following is true, which unfortunately remains only a conjec- 
ture, appearing in several places, most recently in [Aid] and |AS03| . 



Conjecture 1 For every c > and r > 3 the limits 

lim jgfa^ Um E[X(n,r)] 

exist. 

The existence of these limits also implies the convergence to the same limits w.h.p. by applying 
Azuma's inequality, see |.TLR00| for the statement and the applicability of this inequality. 

The limit lim„ EX(n, c)/n is known to exist for c < e and can be computed using the Karp- 
Sipser |KS81j algorithm for maximum matching A4(n, c). We describe the algorithm first in high 
level terms, and then state the result and its implications to maximum independent sets. The 
algorithm proceeds in two stages. In the first stage any leaf v in the graph G(n, c/n) is selected. 
The edge incident to this leaf is selected into a matching, and all the other edges incident to 
the parent of v are deleted. This is repeated until no leaves are left. In the second stage simply 
the largest matching M* in the remaining graph G* C G is selected and added to the matching 
constructed in first stage. It is a simple exercise to prove that the resulting matching is optimal. 
Karp and Sipser prove the following result about the size M. {KS) of the obtained matching. 

Theorem 1 (Karp and Sipser [KS81J) The constructed matching is optimal, A4(KS) =Ai(n,c) 
with probability one. Moreover, 
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The maximum matching satisfies 




1 - 



7*(c) +7**(c) + c 7 *(c)7**(c) 
2 



where 7*(c) zs i/ie smallest solution of the equation x = exp(— cexp(— cx)) and 7**(c) = 
exp(-c 7 *(c)). 

U. When c < e, the equation x = exp(— cexp(— cx)) has a unique solution 7(c), and \G*\ = 
o(n). TTtat is t/ie matching constructed in the first stage is asymptotically optimal as n — > 
00. 

The different behavior for c < e and c > e is called e-cutoff phenomena. Theorem ^ was 
strengthened later by by Aronson, Frieze and Pittel [APF98 , who obtain bounds on convergence 
It is the second part of the theorem above that can be used for the analysis of I(n,c). 
Observe, that in the first stage of the Karp-Sipser algorithm for every selected leaf if one takes 
a parent of this leaf instead of an edge between the leaf and its parent, one obtains a node set 
which is an edge-cover set in the graph G \ G* constructed in the first stage. That is every edge 
in G \ G* is incident to at least one node in this cover. It is a simple exercise to prove that this 
is in fact a minimum node cover in G \ G*. Its complement is a maximum independent set in 
G\G* . Since \G*\ = o(ri) when c < e, one obtains then the following result. 

Corollary 1 When c < e: 



The Karp-Sipser algorithm hinges strongly on working with leaves and thus is not applicable 
to random regular graphs. Moreover, if the edges or the nodes of the graph G(n,c/n) are 
equipped with weights then the Karp-Sipser algorithm clearly can produce a strictly suboptimal 
solution and cannot be used in our setting of weighted nodes and edges. Also, when the edges of 
G r (n) are equipped with weights, the problem of computing maximum weight matching becomes 
non-trivial, as opposed to the unweighted case when the full matching exists w.h.p. 

In a somewhat different domain of extremal combinatorics the following result was established 
by Hopkins and Staton HS82J. A girth of a graph is the size of the smallest cycle. It is shown 
in [HS82 that the size of a largest independent set in an ra-node graph with largest degree 3 and 
large girth is asymptotically at least (7/18)n — o(n). The techniques we employ in this paper 
allow us to improved this lower bound. 



We begin by introducing the key technique for our analysis - recursive distributional equations 
and its fixed point solutions. This technique was introduced by Aldous Ald92 , [AldOTj in 
the context of solving ((2) limit conjecture for random minimal assignment problem, and was 
further developed in Aldous and Steele [AS03j . Steele [Ste02j . Aldous and Bandyopadhyay |AB| . 
Gamarnik |Gam02j . Let W be a non-negative random variable with a distribution function 
F w (t) = F(W < t). We consider four operators T = Tx, n Tr >c , Tm,ti Tm,c acting on the space of 
distribution functions F(t),t > 0, where c > is a fixed constant and r > 2 is a fixed integer. 



(2) 



lim 



E[T(n,c)] _ 27(c) + 7 2 (c) 
n 2 



3 Main results 
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1. Given W distributed according to F w (we write simply W ~ F w ), and given a distribution 
function F = F(t), let Bi, B 2 , . . . , B r ~ F be generated independently. Then T I>r : F — >■ 
F', where F' is the distribution function of B' defined by 

(3) B' = max(0,W- B i)- 

\<i<r 

2. Under the same setting as above, let Bi, . . . , B m ~ F, where m is a random variable 
distributed according to a Poisson distribution with parameter c, independently from W, Bi. 
Then Tj c : F — > F', where F' is the distribution function of B' defined by 

(4) B' = max(0,W- B ^ 

l<i<m 

when m > 1 and B' = W when m = 0. For simplicity we will identify the sum above with 
zero when m = 0. 

3. Let Wi, . . . , W r ~ F w , B 1 ,...,B r ~ F. Then T M>r . : F -> F', where F' is the distribution 
function of F' defined by 

(5) F'= max(0,W / i -F i ). 

l<i<r 



4. Finally, let Wi, . . . ,W m ~ F to , B 1 ,...,B m ~ F, where m ~ Pois(c), independent from 
W{,Bi. Then Tvj c : F — > F', where F' is the distribution function of F' defined by 

(6) F' = max (0, W t - B { ), 

l<i<m 

when m > 1 and F' = when m = 0. Again, for simplicity, we assume that max expression 
above is zero when m = 0. 



A distribution function F is defined to be a fixed point distribution of an operator T if 
T(F) = F. 

We now state the main result of this paper. Recall, that a distribution function F(t) is 
defined to be continuous (atom free) if for every x in its support lim e _ >0 (-^( a; + e ) ~ -^(^ ~ e )) — 0. 
Equivalently, for B ~ F and every x, P(F = x) =0. We use 1{-} to denote the indicator 
function. 

Theorem 2 Let F w be a continuous non-negative distribution function. Forr > 1 if the operator 
Txr-i ^ as a un, iq ue fixed point distribution function F* , then, w.h.p. 

(7) lim X w (n,r) = _ y- Bj>0}]; 

n n 

l<i<r 

where W ~ F„, ; Fj ~ F* ; and W, Fj are independent. When G r (n) is replaced by G(n,c/n), the 
same result holds for T = Tj tC , except the sum in the right-hand side of J7p is Xa<i< m B i an< ^ 
m ~ Pois(c). 
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Finally, the similar results hold for Ai w (n,r) and M. w {n,c) in G r (n) and G(n,c), for T = 
Tm,v-i an d T = Tm,cj respectively, whenever the corresponding operator T is such that T 2 has 
the unique fixed point distribution F* . The corresponding limits are 

(8) lim M^hll = 1 E [ V Wi l{Wi -Bi= max (Wj - BA > 0}], 

n n 2 Kj<r 

l<i<r 

where Wi ~ F w , Bi ~ F* , and 

(9) lim Mw ^ c "> = Ie[V Wi l{Wi -B t = max(iy j - BA > 0}], 

where Wi ~ F w , B { ~ F*, m ~ Pois(c). 

For G = G r (n),G(n,c/n) with the weights on nodes given by a distribution function F w , 
let TJ\f w (n,r), IAf w (n,c) denote the cardinality of the maximum weight independent set in 
G which achieves the maximum weight. In case F w is continuous, the maximum weight in- 
dependent set is uniquely defined, so TJ\f w (n,r),TAf w (n,c) are well-defined as well. Clearly, 
X(n, r) > XJ\f w (n,r),I(n,c) > TJ\f w (n,c). M.N w {n,r) and A4Af w (n,c) are defined similarly. 
When T| r or T| r have the unique fixed point it is also possible to compute asymptotically 
lj\f w (n, r),lj\f w (n, c), MAf w (n, r),MAf w (n, c). 

Corollary 2 Under the setting of Theorem^ w.h.p. 

(10) hm JA ^ (n ' r) Z \{W Y Bi>0}], 

n n — ' 

\<i<T 

where W ~ F w , Bi ~ F* , W, Bi independent, and F* is the unique fixed point distribution of 
T 2 = T| r „ 1; and 

(11) lha MAT w (n,r) = 1 _ > 

where F w , Bi ~ F* , Wi, B^ independent, and F* is the unique fixed point distribution of 

T 2 = Tj^ r _ x and W, Bi are independent, 

Similar results hold for T = Tx, c >Tm,c where again in the sum J2i<j< r we substitute r with 
random m ~ Pois(c). 

The Theorem El is the core result of this paper. It will allow us to obtain several interesting 
corollaries, which we state below. 

Theorem 3 Suppose the weights of the nodes and edges of the graphs G = G r (n) and G = 
G(n,c) are distributed as Exp(l). Then 
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1. Tj r _ 1 has a unique fixed point distribution F* iff r < 4. In this case, w.h.p. 

(12) lim T w (n,r) = (l-b)(r~rb + 2b + 2) 

n n 4 

where b is the unique solution ofb=l — (^^) r_1 . In particular, w.h.p. 

(13) lim J " (n ' 2) =g lim 6077, lim 5^ « .4974, 

n n 3 « n n n 

£. T| c has a unique fixed point distribution F* iff c <2e. In this case, w.h.p. 

(14) lim X>,c)^ c(l-b) 

n n 4 

where b is the unique solution of 1 — b = e - ^ 1-6 ). /n particular, when c = 2e, this limit is 
ss .5517. 

5. r-1 /ias a unique fixed point F* for every r > 2. Moreover, w.h.p. 
(15) 

lim wK ' ; =r(6 r - 1 + l) / te-*(l-e-*(l-6))) r - 1 c/t-r / - e -*(l -6))) 2r " 2 rft, 

" io Jo 

where b is the unique solution ofb — 1— r ]j^ ■ 
4- Tjtf c a unique fixed point F* for all c> 0. Moreover, for every c > w.h.p. 

(16) lim M^c) = C (eCb _ c + 1) r te -t-c(l- b )e-t dt _C /•°° te _ i _ 2c(1 _ b)e - +c(1 _ fe)2e - 2t ^ 

n n 2 y 2 J 

where b is the unique solution of 1 — e~ c6 = c(l — b) 2 . 

The result above generalizes to the case when Exp(l) is replaced by Exp(/x) for any fi > 0, 
since for any a > and W ~ Exp(l), al^ ~ Exp(l/a). The expression in (|15|) involving 
integrals is similar to the one found in |AS03| for maximum weight matching on a tree. It is a 
pleasant surprise, though, that the answers for independent sets are derived in closed form. 

Part 2 of Theorem El leads to an interesting phase transition behavior. For F w = Exp(l) 
our result says that the value of c = 2e is a phase transition point for the operator Tj c where 
for c < 2e the operator has a unique fixed point distribution, but for c > 2e the fixed point 
distribution is not unique. Contrast this with e-cutoff phenomena described above. It turns 
out (see Theorems EH below) that this phase transition is directly related to some long range 
independence/dependence property of the maximum weight independent sets in the underlying 
graph G(n,c/n). Curiously, no such phase transition occurs for maximum weight matchings. 

As a sanity check, let us show directly the validity of ()12j) for the case r = 1. The answer 
given by the formula is 3/4. It is, though, easy to compute I w (n, 1) exactly. The graph is a 
collection of n/2 isolated edges. For each edge we simply select the the incident node with larger 
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weight, which by memory less property of Exp is 1/2 + 1 = 3/2. The limit T w {n, l)/n — > 3/4 then 
checks. 

Note that G2(n) can be represented as a collection of disjoint cycles on n elements. W.h.p. 
these cycles have decreasing lengths, starting from 0(77.2). Theorems E] and |3] also hold also 
for a simpler model of a 2-regular graph - n-cycle. The n nodes 0,1,..., n — 1 and edges 
(0, 1), . . . , (n — 2, n — 1), (n — 1, 0) of this cycle are assumed to have weights distributed according 
to some distribution function F w . Let I w (n, cycle) and M w (n, cycle) denote respectively the 
maximum weight of an independent set and the maximum weight of a matching. 

Corollary 3 Suppose T| 2 has a unique fixed point. Then |7p ; hold when r = 2 and 
X w (n, cycle) replaces I w (n, 2) and A4 w (n, cycle) replaces M. w {n, 2). When F w = Exp(l), w.h.p. 

(17) lim X w (n, cycle) _ ^ M w (n, cycle) _ 2 

n n n n 3 

We can also compute the cardinality of independent sets which achieve the maximum weight, 
when F w = Exp(l). 

Corollary 4 When F w = Exp(l), w.h.p. 

( 18 ) , im inf « > „ m ZV. (n ,2) = 4 Hm inf XK3) > hm ZV.( B ,3) M 3923 _ 

nn«n 9 nn«n 

l iminf ^j)>l, m ^("- 4 ) K .3533. 

n n n n 

The last part of the corollary above implies a new lower bound on the size of the largest inde- 
pendent sets in a random 4-regular graph. 

Corollary 5 The cardinality of the largest independent set in G±(n) is at least .3533n — o(n). 

We note, however, that for r = 3 our lower bound .3923n — o(n) is weaker than A32n — o(n) 
established in TW\}. Yet, it is still very useful since the approach allows us the improve the 
following result of Hopkins and Staton |HS82j . Let Q (n, r, d) denote the class of all (non-random) 
graphs on n nodes, with maximum degree r and girth at least d. For any G G Q(n,r,d) let 
T{G) denote the size of the largest independent set in G. Hopkins and Staton proved that 
liminf n) dmin Gg g( ni 3 irf )Z(G)/n > 7/18 ~ .3887. Our techniques allow us to obtain the following 
improvement. 

Theorem 4 

lim inf min ?2f2 > .3923. 

n,d G&g(n,3,d) 71 

An important implication of the uniqueness of the fixed point distribution of T 2 for the types 
of T described above, is that the uniqueness implies a certain long-range independence properties 
of the structures we consider. The following theorem makes this notion precise. While the 
theorem is not used directly in this paper, we believe it is interesting by itself. Below G is again 
one of the graphs G r (n) or G(n, c). Let E denote the (random) edge set of G. 
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Theorem 5 Let T be one of the four operators J3j), J3J) ; (QJ) witt respect to some continuous 
distribution function F w , and let C w {n) = T w (n,r),T w (n,c),J\A w (n,r) or Ai w (n,c). Denote by 
0{n) the subset of [n] or E (depending on a context), which achieves C w {n). Select two elements 
i,j of [n] or E uniformly at random. IfT 2 has a unique fixed point distribution, then 

(19) G 0(n)) -> F(i G 0(n))F(j e 0(n)), 

as n — > oo. 

Recall that for each of the four objects \C{n)\ = 0(n) w.h.p. As a result the values P(t> G 
0(n)) do not vanish and the first part of the theorem above does have a non-trivial content. In 
fact we will show a much stronger result stating that, when T 2 has a unique fixed point, the 
event v G 0(n) is almost independent of the entire graph G outside a depth-d graph-theoretic 
neighborhood of v, when d is a sufficiently large constant integer. 

4 Fixed points of the operator T 2 and the long-range in- 
dependence 

4.1 Maximum weight independent sets and matchings in trees. Fixed 
points of T 2 and the bonus function 

We start by analyzing operator T - one of the four operators introduced in the previous section. 
Given two distribution functions Fx, F 2 defined on [0, oo), we say that F 2 stochastically dominates 
F\ and write F\ -< F 2 if F\(t) > F 2 (t) for every t > 0. A sequence of distribution functions F n is 
defined to converge weakly to a distribution function F (written F n =>- F) if lim n F n (t) = F(t) 
for every t which is a point of continuity of F. 

Lemma 6 The operators T = Tx r , Tj jC , Tj^ >r , Tj c are continuous with respect to the weak con- 
vergence. That is, given a sequence of distributions F, F s , s = 1, 2, . . . , if F s F then T(F S ) =>• 
T(F). 

Proof : The proof is almost immediate for T = Tj r . Since the summation, subtraction and 
max operations are continuous functions then the assertion holds. Here we use the fact that for 
any continuous function f,F n ^F implies f{F n ) =^ f{F) (see Continuous Mapping Theorem in 
|Dur96j ). The proof for the case T = Tj jC is slightly more subtle since we are dealing with sum of 
randomly many elements. Let x[ , X% , . . . , Xm' be distributed independently according to F n 
and let m be random variable with Pois(c) distribution. Let X\, . . . , X m be distributed according 
to F. Define X^ = max(0, W - Ei<*< m x< t ] ) and X = max(0, W - Ei<;< m X r)- Fix arbitrary 
e > and m such that m does not exceed m with probability at least 1 — e. For any fixed t > 

|P(X (n) <t)- P(X (n) < t\m < m )P(m < m )| < e, 

and 

\F{X <t)- F(X < t\m < m )P(m < m )\ < e, 
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Now, we have |P(XW < t\ th < "T,o ) — P(X < t\m < mo) | < e for sufficiently large n since F n F 
and m is conditioned to be at most m . Combining, we obtain |P(X^ < t) — F(X <t)\ < 3e for 
sufficiently large n. This completes the proof for T = Tj c . The proofs for T = Tm,t, T = Tm, c 
are similar. □ 
Let denote (for simplicity) the distribution function of a random variable X which is zero 
w.p.f. Let W r = max!<j< r Wi where Wi ~ F w are independent. Let also W c = max!<j< m Wi, 
where Wi ~ F w are independent and m ~ Pois(c). Denote by F w>r and F w>c the distribution 
functions of W r and W c , respectively. 

Proposition 1 Fix r > l,c > ; a distribution function F w (t),t > and T = Tj r or Tj c . ^4s 
s — > oo, t/ie two sequences of distributions T 2s (0),T 2s (F w ) weakly converge respectively to some 
distribution functions F**,F** which are fixed points of the operator T 2 . For any distribution 
function F = F (t),t > 0, T 2s (0) -< T 2s (F ) -< T 2s (F w ) and T 2s+1 (F W ) -< T 2s+1 (F ) -< T 2s+1 (0) 
for all s = 1, 2, . . .. 

If the operator T is such that T 2 has a unique fixed point (and F** = F** = F*), then for 
any distribution function F = F (t),t > 0, T S (F ), s = 1,2,... converges to F* as s —>■ oo. In 
particular, T s (0),T s (F w ) — > F* . Moreover, F* is also the unique fixed point ofT. 

When T = T M , r or T = T MjC the same result holds with F w ^ r and F w ^ c respectively replacing 

F 

Proof : Let F = F(t),t > be any distribution function corresponding to any non- negative 
random variable. It follows immediately from the definitions that 

(20) Tx )T {F) -< F w , Ti jC (F) -< F w , Tm,t{F) -< F w>r , Tm,c{F) -< F w>c . 

Observe that the each of the four operators T above is anti-monotone. That is if F\ -< F2 
for some distribution functions F±,F2, then T(Fi) >- T(i^). Applying this twice we obtain 
T 2 (Fi) -< T 2 (F 2 ). Then -< T 2 (0) -< T 4 (0) -<•••, and this sequence weakly converges to some 
function F^ = F(t)^ < 1, since the values of T 2s (0) at each fixed t are decreasing and are 
bounded below by 0. Note that F^(t) is a non-decreasing function of t, since this is the case for 
each distribution T 2s (0). Finally, from (J2TJ|) we have that lirn^oo F**(t) = 1. Thus F^(t) is a 
distribution function. As T 2s (0) =^ F** and by LemmaEl T is continuous, we have T 2 (F„) = F**. 

We now fix T = Tj jr . The proof of the other three cases is very similar. From (J2(J|) . by taking 
F = Tj t r(F w ) we obtain T 2 (F W ) -< F w . Then, by monotonicity of T 2 we obtain F w y T 2 (F W ) y 
• • ■ y T 2s (F w ) and the sequence T 2s (F w ) converges weakly to some function F**. We repeat the 
arguments above to show that F** is actually a distribution function. Again applying Lemma 
we conclude T 2 (F**) = F** . 

Suppose now T 2 has a unique fixed point F* = F** = F**. For any distribution function 
F we have -< T(F ),T 2 (F ) -< F w , where again we use (|2*U|) . Applying the monotonicity 
T 2s (0) -< T 2s+1 (F ),T 2s+2 (F ) -< T 2s (F w ). But since T 2s (0), T 2s (F w ) F*, then T S (F ) F*. 
In particular, by taking F = and F = F w , we obtain T s (0), T S (F W ) — > F*. Finally, we obtain 
T 2s (T(F*)) -> F*. But T 2s (T(F*)) = T(T 2s (F*)) = T(F*). Thus T(F*) = F* . Clearly it is the 
unique fixed point of T since any fixed point of T is also a fixed point of T 2 . □ 

We now switch to analyzing the maximum weight independent set problem on a tree. The 
derivation here repeats the development in |ASQ3j for maximum weight matching in random 



trees. We highlight important differences where appropriate. 
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Suppose we have a (non-random) finite tree H with nodes 0,1, . . . , h = \H\ — 1, with a fixed 
root 0. The nodes of this tree are equipped with some (non-random) weights W , W\, . . . , Wh > 0. 
For any node i G H, let H{i) denote the subtree rooted at i consisting of all the descendants 
of i. In particular, H(0) = H. Let 1h{%) denote the maximum weight of an independent set in 
H(i) and let B H ^ = Tu(i) — J2j^H(j), where the sum runs over nodes j which are children of 
i. If i has no children then Bh(%) is simply Th(i) — W*. Observe, that Bhu) is also a difference 
between Th(%) an d the maximum weight of an independent set in H(i), which is not allowed to 
use node i. Clearly, < B H ^ < Wi. The value B H ^ was considered in |AS03j in the context 
of maximum weight matchings and was referred to as a bonus of a node % in tree H{i). W.l.g. 
denote by 1, . . . ,m the children of the root node 0. 

Lemma 7 

(21) B m = max(0, W - ^ B H{i) ). 

l<i<m 

Moreover, if W > Ylii<i<m B H{i) (that is if B h iq\ > 0) then the maximum weight independent 
set must contain node 0. If Wo < J2i<i<m Bu(i) then the maximum weight independent set does 
not contain the node 0. 

Remark : There might be several independent sets in H which achieve 1h- The second part 
of the lemma refers to any independent set achieving maximum weight. Also, the statement of 
the lemma applies as well to every node i with respect to its tree H(i). That is let ji, ■ ■ ■ ,ji be 
the children of a node i. Then the lemma claims B H t^ = max(0, W — J2i<i<i ^H(ji))- 

Proof : Consider an independent set V C H which achieves the maximum weight. We 
take an arbitrary such in case there are many. If G V then i ^ V for 1 < i < m. Then V is 
obtained by taking maximum independent sets Vi in H(i) such that i ^ V(i). By definition the 
weight of such Vi is lH(i) — B H ^y On the other hand, if ^ V, then a maximum weight of an 
independent set in H is obtained simply as ^2 l<i<m ^H(i)- We conclude that 

(22) 1 H = max(^ + ^ ( J "« ~ B n®)> Ih ^' 

l<i<m l<i<m 

Recall that B H = X H — J2i<i< m -^H(i)- Subtracting ^ 1<j<m %) from both sides of (|2^|) we 
obtain (|21jl . The second part of the lemma follows directly from the discussion above. □ 
A similar development is possible for maximum weight matching. Suppose the edges of the 
tree H are equipped with weights Wij. Let M.H{i) denote the maximum weight of a matching in 
H(i), and let Bh(i) denote the difference between M.h{i) and the maximum weight of a matching 
in H(i) which is not allowed to include any edge in H(i) incident to i. Again 1,2, ... ,m are 
assumed to be the children of the root 0. 

Lemma 8 

(23) B H = max(0, max (W 0>i - B H[i) )). 

l<i<m 

Moreover, ifWo t i — B H ^ > Wq^ —B h ^/) for alii' ^ i and Wo,i — B H ^ > 0, then every maximum 
weight matching contains edge (0, i). IfWo^ — Bu^ < for alii = \, ... ,m, then every maximum 
weight matching does not contain any edge incident to 0. 
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Proof : The proof is very similar to the one of Lemma [7| If a maximum weight matching 
contains edge (0, i) then its weight is Wq^ plus maximum weight of a matching in H(i) with node 
i excluded, plus the sum of the maximum weights of matchings in Tn(j), 1 < j < m,j ^ i. If a 
maximum weight matching contains none of the edges (0, i), 1 < i < m, then its weight is simply 
sum of maximum weight of matchings in H{i) for all i — 1, . . . , m. We obtain 

(24) M H = max( max (W 0)i + M H {i) - B H{i) + V M H (j)), Y] M H (j))- 

j^i l<i<m 

Subtracting Y^i<j< m ^ H {i) fr° m both sides of (|2^j) we obtain (|2^j) . The proof of the second part 
follows immediately from the discussion above. □ 

4.2 Long-range independence 

We now consider trees H of specific types. Given integers r > 3, d > 2 let H r (d) denote an 
r-regular finite tree with depth d. The root node has degree r — 1, all the nodes at distance 
> 1, < d — 1 from the root have outdegree r — 1, and all the nodes at distance d from are 
leaves. (Usually, in the definition of an r-regular tree, the root node is assumed to have degree 
r, not r — 1. The slight distinction here is done for convenience.) Also, given a constant c > 0, a 
Poisson tree H(c, d) with parameter c and depth d is constructed as follows. The root node has a 
degree which is a random variable distributed according to Pois(c) distribution. All the children 
of have outdegrees which are also random, distributed according to Pois(c). In particular, the 
children of have total degrees l + Pois(c). Similarly, children of children of also have outdegree 
Pois(c), etc. We continue this process until either the process stops at some depth d! < d, where 
no nodes in level d' have any children, or until we reach level d. In this case all the children of 
the nodes in level d are deleted and the nodes in level d become leaves. We obtain a tree with 
depth < d. We call this a depth-rf Poisson tree. 

Let H = H r (d) or H (c, d). Suppose the nodes and the edges of H are equipped with weights 
Wi,Wij, which are generated at random independently using a distribution function F w . Fix 
any infinite sequences w = (wi, w 2 , ■ ■ .) G [0, oo)°° and b = (bi,b 2 ,-.-) G {0,1}°°. For every 
i — 1, 2, . . . , d let il, i2, . . . , iji denote the nodes of H in level % (if any exist for H(c, d)). When 
H = H r (d),ji = (r — 1)% of course. Let (I\{b, w)) denote the maximum weight of an independent 
set V in H such that the nodes dj with bj = 1 are conditioned to be in V, nodes dj with bj = 
are conditioned not to be in V, and the weights of nodes dj are conditioned to be equal to Wj 
for j = 1, . . . ,jd- That is we are looking for maximum weight of an independent set among 
those which contain depth d leaves with bj = 1, do not contain depth d leaves with bj = 0, and 
with the weights of the leaves deterministically set by w. For brevity we call it the maximum 
weight of an independent set with boundary condition (b, w). For the case H = H(c,d), the 
boundary condition is simply absent when the tree does not contain any nodes in the last level d. 
(TH(ij)\{b,w)) are defined similarly for the subtrees H(ij) spanned by nodes ij in level i: given 
again b, w, let (A4\(b, w)) and (A4H(ij)\(b, w)) denote, respectively, the maximum weight of a 
matching E in H and Hij, such that the edges incident to nodes dj are conditioned to be in E 
when bj = 1, edges incident to nodes dj are conditioned not to be in D when with bj = 0, and 
the weights of the edges incident to nodes dj are conditioned to be equal to Wj,j = 1, . . . ,jd (of 
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course, we refer to edges between nodes in levels d — 1 and d as there is only one edge per each 
node in level d). 

For the case of independent sets, let (B\(b, w)) denote the bonus of the root node given the 
boundary condition (b, w). Namely, 

(B\(b,w)) = (I\(b,w))- 

i<j<h 

For the case of matchings, let (B\(jb,w)) also denote the bonus of the root node given the 
boundary condition (b, w). Namely, 

(B\(b,w)) = (M\(b,w))- (M H{l3) \(b,w)). 

i<j<h 

It should be always clear from the context whether B is taken with respect to independent sets 
or matchings. 

The following theorem establishes the crucial long-range independence property for the maxi- 
mum weight independent sets and matchings in trees H = H r (d), H(c, d) when the corresponding 
operator T 2 has a unique fixed point. It establishes that the distribution of the bonus (B\(b,w)) 
of the root asymptotically is independent of the boundary condition (b, w) as d becomes large. 
Recall our convention that denotes the distribution of a random variable which is zero to one 
with probability one. 

Theorem 9 Given a distribution function F w and a regular tree H = H r (d) let T = T%^. Then 
for every t > and b, w 

(25) T d -\0)(t) < P((B\(b,w)) <t)< T d -\F w ){t) 
when d is odd, and 

(26) T d -\F w )(t) < n(B\(b,w)) <t)< T d -\0)(t) 

when d is even. Suppose in addition the operator T 2 = Tj r _ 1 has the unique fixed distribution 
F*. Then 



(27) sup ¥((B\(b,w)) < t)-F*(t) 

b,w 



0. 



as d —>■ oo. Similar assertion holds for T = T^ r and for H = H(c,d) with T = Tj ;C and 
T = Tm, c - For the cases T = Tm,t and Tm, c , F w is replaced with F WjT and F WtC respectively. 

Before we prove the proposition above, which is essentially the key result of the paper, let 
us compare it with the developments in ([AS03 ). In that paper maximum weight matching is 
considered on an n-node tree, drawn independently and uniformly from the space of all n n ~ 2 
labelled trees. The notion of a bonus is introduced and the recursion is derived. However, 
since a tree structure is assumed to begin with, there is no need to consider the boundary 
conditions (6, w) . Here we avoid the difficulty of the non-tree structure by proving the long- 
range independence property via the uniqueness of fixed points of T 2 . 



14 



Proof : We prove ([250 and (|26|) for independent sets in H = H r (d). The proof of (|27|1 will 
then be almost immediate. The proofs for other cases is similar. We will indicate the differences 
where appropriate. The proof of (}25|) proceeds by induction in d. Suppose d — 1. Then we have 



(B\(b,w)) = max(0, W ( 



i:l<i<r—l,bu=l 



Then (F|(6, u>)) -< Fu,. Trivially -< (B\(b,w)). This establishes the bound for the case d = 1. 
Suppose the assumption holds for d — 1 and d is odd. We have 



By the inductive assumption the distribution of each of (Bu\(b,w)) dominates T d ~ 2 (F w ) and is 
dominated by T d ~ 2 (0). Applying this to (|28|) we obtain (|25|) . The case d is even is considered 
similarly. This completes the induction and establishes (|25|) and (|26p . Now if T 2 has the unique 
fixed point distribution F* then by Proposition ^ T d (0), T d (F w ) converge to F* as d — ► oo, and 
we obtain (|2*7j) . The proof for the case of Poisson tree H(c,d) instead of H r (d) is pretty much 
identical. For the proof of the same result for maximum weight matching we just use F Wjr and 
F w c instead of F w , just as we did in the proof of Proposition [T] □ 
While it is not important for the further results in this paper, it is interesting that the 
uniqueness of the solution T 2 (F) = F is the tight condition for (J2*7)l . as the following theorem 
indicates. 

Theorem 10 Suppose the operator T 2 has more than one fixed point distributions F* . Then for 
every such F* 



Proof : As usual we start with T = Tj r . The proofs for other cases are similar, we highlight 
the differences where appropriate. Let F„ and F** be distributions introduced in Proposition 
[3 The non-uniqueness of the fixed point of T 2 implies using Proposition Q that F** ^ F**. For 
every node j in layer d (last layer) of H r (d) set bj = 0,Wj = 0. In particular, the bonus Bj of 
each such node is zero. Then for every node in layer d — 1 its bonus is given by the distribution 
T(0) = F w , the bonus of each node in layer d — 2 has distribution T 2 (0), etc. The root node has 
bonus with distribution T d (0). When d is an even number diverging to infinity, from Proposition 
HJ T d (0) converges weakly to F**. Thus the distribution of (B\(b,w)) converges to F**. If F*, 
a fixed point of T 2 , is distinct from F**, then we obtain that (f2HJ) holds. Suppose, on the other 
hand, F* = F**. We claim that T(F^) ^ F**. Assuming this is the case we consider the same 
boundary condition b = w = but take d to be odd integer diverging to infinity. Then the bonus 
of the root converges in distribution to T(F^) ^ F** and (}2T?|) is shown again. 

Assume T(F^) = F**. Recall from the first part of Proposition ^ that for every distribution 
F , T 2s+1 (F ) -< T 2s+1 (0). Taking F = T(F**) we obtain F** = T 2s+2 (F**) -< T 2s+1 (0). Taking 
s — > oo, F** -< T(F^) = F«. But from Proposition [T] F** -< F**, and, as a result F** = F** 
implying (again using Proposition ^) T 2 has the unique fixed point. We obtained a contradiction. 



(28) 




i:l<i<r— 1 



(29) 



liminfsup F((B\(b, w)) < t) - F*(t) > 0. 



b,w 



□ 
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5 Applications to maximum weight independent sets and 
matchings in G r (n) and G(n,c/n) 

5.1 Long-range independence in G r (n),G(n,c/n) 

The goal of the current section is to demonstrate that Theorem |U] allows us to reduce the com- 
putation of the maximum weight independent set and the maximum weight matching in random 
graphs to a much simpler problem of finding those in trees. We highlight this key message of 
the paper as the following local optimality property: if the operator T 2 corresponding to a max- 
imum weight combinatorial object (independent set or matching) in a sparse random graph has 
a unique fixed point, then for a randomly selected node (edge) of the graph, the event "the node 
(edge) belongs to the optimal object" and the distribution of the node (edge) weight, conditioned 
that it does, asymptotically depends only on the constant size neighborhood of the node and is 
independent from the rest of the graph. In other words, when T 2 has a unique fixed point, the 
maximum weight independent sets and matchings exhibit a long-range independence property. 

Our hope is that similar local optimality can be established for other random combinatorial 
structures. 

Proof : [Proof of Theorem Ej. Again we start by proving the result for 1 w (n,r). The 
proofs for other three objects is similar and we highlight some differences in the end. Let 
V r C G r (n) denote the independent set which achieves the maximum weight X w (n,r). It is 
unique by continuity of the distribution F w . Consider a randomly selected node of the graph G, 
which, w.l.g., we may assume is node 0. By symmetry we have E[Z M (n,r)] = riE[Wol{0 € V r }}. 
Let us fix a large positive integer d, which is a constant independent from n, and let H (d) denote 
the depth-d neighborhood of 0. That is H is the collection of nodes in G which are connected 
to by paths with length < d. It is well known that, w.h.p. as n — > oo, H is a depth-c? r- 
regular tree, |.ILR00j . except in this case, unlike in Subsection \A.2\ the root node has outdegree 
r and the remaining non-leaf nodes have outdegree r — 1. Let dH denote the leaves of this 
tree (level d). Fix any binary vector b with dimension \dH\ + \G\H\, and any non-negative 
vector w also with dimension \dH\ + \G\H\. Assume the vector b is such that if two nodes in 
dH U (G \ H) are connected by an edge, only one of these two nodes can have the corresponding 
component of b equal to 1. That is, the nodes marked 1 by b correspond to some independent set 
in dH U (G \ H) . Consider the problem of finding the maximum weight independent set V r in G 
when the weights of the nodes in dH U(G\ H) are conditioned to be w, nodes i G dH U (G \ H) 
with the corresponding component of b equal to one are conditioned to belong to V r and the 
remaining nodes in dH U (G\H) are conditioned not to belong to V r . Notation- wise, we consider 
the value (I w (n,r)\(b,w)). In particular, we need to select to maximum weight independent set 
in the tree H, which is consistent with conditioning (b,w). Naturally, the consistency needs to 
be checked only across the boundary dH. 

Let (B\(b, w)) and (Bi\(b,w)), 1 < i < r denote the bonus of the node and the bonuses of 
its neighboring nodes 1,2, ... ,r, respectively. By Lemma [7[ we have [B\(b, w)) = max(0, Wo — 
^2 1<i<r (Bi\(b, w))). Since the distribution function F w of W is continuous, then F* is continuous 
as well and, as a result, W — J2i<i< r (Bi\(^w)) = with probability zero. Then applying the 
second part of Lemma node belongs to the maximum weight independent set if and only if 
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B = Wo - Ei<i< r (^l(M)) > 0. Therefore 

E[J ^' r)] = £E[W i{Wb - (s*|(5, «)) > 0}]P(M). 

6,ui l<i<r 

But by Theorem ITUl for every &, w the distribution of (-Bj|(5, w)) converges to F*, the unique 
fixed point of T| r _ 1 as d becomes large. We conclude 

ljm E^(n,r)] = E[wW _ ^ 5i > o}], 

n 77, ^ — ' 

l<i<r 

where ~ F*. This completes the proof for the maximum weight independent set in G r (n). 
When the graph G(n,c/n) is considered the proof is very similar, we just use the fact that H 
- depth d neighborhood of approaches in distribution a Poisson tree |Spe94| . The proofs for 
maximum weight matchings are similar, we use Lemma |S1 instead of Lemma [7| □ 

Proof : [Proof of Corollary |2] In fact we have proved this result already en route of proving 
Theorem El above. We have shown that given b, w and conditioned on H being a tree a fixed node 
belongs to the maximum weight independent set iff (B\(b,w)) = Wo — J2i<i< r (Bi\(b,w)) > 0. 
Repeating the proof of Theorem f%IAf w (n, r)/n = E[1{0 G V r }] -> E[1{W - Ei<i<r B ^ where 
Bi ~ F*. □ 

Proof : [Proof of Theorem^ The essential ingredients for the proof of this results are already 
established in the proof of Theorem El above. We start with X w (n,r). The proofs for other cases 
are very similar and we omit them. Let O = 0(n,r) C [n] be the independent set achieving 
the maximum weight. Fix two nodes i,j G [n] and arbitrary e > 0. Let H = H(i,d) denote 
the depth-<i graph-theoretic neighborhood of % in the graph G = G r (n) and let dH denote the 
boundary of H - the nodes of H at distance d from i. Let St denote the event that if is a 
(r-regular) tree. From the theory of random regular graphs |JLR00j . P(£t) - ► 1 as n — > oo. As 
a result ¥(St) > 1 — e for all n > no for some no = no(d) (note the dependence on d). Fix any 
realization of (G \ H) U dH together with the realization of the weights w of nodes in G \ H U OH 
and indicators b of whether the nodes belong to O. As far as deciding which nodes of H\dH are 
in O and in particular whether node i belongs to O only the restriction of w, b to dH is relevant. 
Conditioning on the event St denote by (B\(b,w)) the bonus of i in H (for completeness define 
(B\(b,w)) to be zero when the event St does not hold). Applying Lemma [7| and the continuity 
of F w , we have i G O iff (B\(b,w)) > 0. Applying TheoremU |P((S|(6, w)) > 0) -¥(B > 0)| < e 
for all d > d (e) for some d (e), and for any w, b, where B ~ F* and F* is the unique fixed point 
of T 2 = Tl r _ v Thus 

|P(z G O) -¥(B > 0)| < ^|P(i eO|£^,G,tw,6)-P(S> 0)|P(G,tw,6) 

+ |P(z G 0\S T ) - f(B > Q)\¥(Sr) 
(30) < 2e, 

whenever d > do(e) and n > no(d), where G denotes generically a realization of the subgraph 
(G\H)U dH. Observe that since \H\ < 1 + r + (r - l) 2 H h (r - l) d then P(j <£H)>l-e 
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for all n > n\{d) for some n\(d). Then 

|P(z,j G 0) -P(i G 0)P(j G 0)| 

(31) < |P(z G 0|£ T ,j ^,J6 0)P(£ T , j i H,j G 0) -P(z G 0)P(j G 0)| 

(32) + P(i G 0,£ T ,j G HHO) G 0,£ T ,j G 0) 

The event j ^ H, j E O is completely described by the realizations G, w, b. Since |P((i?|(6, w)) > 
0) - P(£ > 0)| < e for all d > d (e) then 

(33) |P(z G 0|£ T , j <£H,jeO)- P(S > 0)| < e, 

for all d > d (e). Also P(j G 0) = P(£r,j £ H,j e O) + ¥{£ T± j E H f] O) + F(£r,j G 0). 
But P(£ T , j 6ffn0)<P(jGF)<e for all n > m(d) and P(£r, j G 0) < F(£ T ) < e for all 
n > n (d(e)). As a result 

(34) |P(j G 0) -P(5 T ,j £ F,j G 0)| < 2e 
whenever n > max(n (c?(e)), ni(c?(e))). Combining (jB^j) . (jSHj) and (jMj) we obtain 

P(z G 0|£ T ,j ^ H,j E 0)F(£ T ,j <£H,jE 0) -P(z G 0)P(j G 0) 

< (P(i G 0) + 3e)(P(j G 0) + 2e) - P(i G 0)P(j G 0) 

< 5e + e 2 < 6e 

Similarly we show 

P(i G 0|£ T , j ^,J£ 0)V(£r,j <£HJ G 0) - P(i G 0)P(j G 0) > -6e. 

We conclude that the value in (JBTj) is bounded by 6e. Each summand in (jB^j) is bounded by e 
since P(j G H) when n > rii(d(e)) and P(£t) < e when n > n (d(e)). We conclude that whenever 
n > max(n (rf(e)),m(d(e))), |P(z, j G 0) -P(z G 0)P(j G 0)| < 8e. This completes the proof of 
the theorem. □ 



5.2 Computation of limits. Exponentially distributed weights 

We prove Theorem El in this subsection. Thanks to Theorem we can focus on proving the 
uniqueness and computing the fixed points of the operator T 2 . As usual we start with max- 
imum weight independent set in G r (n). The analysis of other cases is similar and will follow 
immediately. The calculations are similar to the ones in [AS03 performed for maximum weight 
matching in random trees. The difference is that we have to compute the fixed point of the 
operator T 2 and not just T. 

Proof : [Proof of Theorem |3J 

• Independent sets in G r (n). Let F* denote any fixed point distribution of T 2 = Tj r _ 1 (at 
least one exists by Proposition [TJ), and let B ~ F*. Then B = max(W — J2i<i< r -i ^*)> ^ ~ 
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Exp(l),S< ~ T(F*). Similarly, if B ~ T(F*), then B = max(W - Ei<i< r -i B i)^ W ~ 
Exp(l), Bi~ F*. Then for any t > 

P(£>t)=P(W> B< + *) = e _t P(W > ^ A), 

l<i<r-l l<i<r-l 

and similarly 

P(S >t) = e-^iW > B ^ 

\<i<r-l 

where we use the memoryless property of the exponential distribution. Let b = F(B = 
0),b = P(J3 = 0) where B ~ F* and I? ~ T(F*). Our next goal is computing 6 and b. 
From above we obtain 

(35) F(B > t) = e _ *(l - 6), P(S > t) = e -t (l - 6), 
implying 

(36) P(S > t|£ > 0) ~ Exp(l), P(B > t\B > 0) ~ Exp(l). 



Then 

(37) 



e -t P( Bi - l ) dt 

i i i— i Ki<r-1 



In order to compute P(Xa<t<r-i ^ — we condition on j < r — 1 terms 5j out of r — 1 
being equal to zero, and the rest positive. This occurs with probability (^^(l - S) r - 1_J '. 

When j < r — 1, the sum of r — 1 — j non-zero terms 5j has an Erlang distribution with 
parameter r — 1 — j (sum of r — 1 — j independent random variables distributed as Exp(l)). 
The density function of this distribution is f(z) = (f_ 2 _jy e ~ z an d the tail probability is 

roc 7 r—2—j j-i 

(38) P( ■ >t)= -e~ z dz = V -e- % 

Jt (r-2-j)\ ^ i 



and 



poo 1 1 

/ ■ >t)dt= y: ^ = 1 -^r- 



0<i<r-2-j 

Combining with (|S7|) and interchanging integration and summation, we obtain 



(39) t= E f r /W-«r 1 - J (i-2^ 3 ) = i-( 1 7^) 



0<?'<r-l 



r-1 



Similar calculations lead to 6 = 1 — (^j^) r 1 . Combining 



1 + ^x r _l 



(40) 6 = /(/(&)), where f(x) = 1 - 
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Lemma 11 The equation MUj) has a unique solution b* within the range b G [0, 1] iffr < 4. 
In this case b* is also the unique solution of f(b) = b. 

Figures 1 and 2 below show the graphs of /(/(&)) for the cases r = 4 and r = 8. The first 
corresponds to the case of the unique solution. In the second case there are more than one 
solution. 

Proof : Note that for every r > 1 the equation x = f(x) has exactly one solution in [0, 1] 
since / is a strictly decreasing function and /(0) = 1 — l/2 r_1 > /(l) = 0. This solution 
is also a solution to x — /(/(#)). 

We now prove the uniqueness for 2 < r < 4 and non-uniqueness for r > 4. Let r < 4. We 
claim that for all x £ [0,1]. 

(41) df(f(x)) = ( r -iy^ + f(x) y-^l + xy-z < ^ 

This would imply that f(f(x)) — x is a strictly decreasing function and is equal to zero in 
at most one point. (J4T]) is equivalent to 

2r-2 

(42) (l + f(x))(l+x)< - 5 . 



Let us check the validity of this inequality at the end points x — 0, 1. Since /(0) < 1, /(l) = 
0, the left hand side in both points is at most 2 which is strictly smaller than the right-hand 
side for r = 2, 3, 4, as it is easily checked. It remains to check the inequality at the points 
where the derivative of the function g(x) = (1 + f(x))(l + x) vanishes. We have 



r — 1 (\ + x\ r - 2 /l + xy- 1 r — l/l+x\ r - 2 

1 — ( — I +1- — — -x— 



g(x) = 1 + f{x) + f(x) + xf(x 

— 1/1 + £ V'- 2 

— (— ) + 

/1 + xy- 1 (1 + x)(r — 1) fl + x\ r 

\~ 

,'l + x\ r - 1 

2 — r 



Thus g(x) = in exactly one point x = 2(-)'- 1 — 1. We check that holds in this point 
for r = 2, 3, 4. 

We now prove the non-uniqueness of the solution to (|40|) when r > 5. Let 6* denote the 
unique fixed point of /(&*) = 6*. Clearly /(/(&*)) = We claim that 

(43) 



> 1. 



This implies the result since we get that for e sufficiently small, f(f(b* — e)) < b* — e. But 
/(/(0)) > 0. Therefore there exists a different fixed point of ()40|) in the interval (0, b*). To 
show (14*31 note that 



df(f(x)) (r - l) 2 ^1 + f(b*)y-2{ 1 + 6* y- 2 _ (r - l) 2 f 1 + b*\^-4 

dx x=b* 4 
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Figure 1: /(/(&)) = b has one solution when r = 4. 



We need to show that 

(44) (r-l^l+^V^ > ^ 



which is equivalent to b* > 2(^j) r - 2 — 1 = b(r). We claim that b(r) < f(b(r)) = 1 — 
( 1+ ^ r - > ) r ~ 1 . Since 6* = /(&*) and / is a decreasing function, this would imply b{r) < b* or 
dSJ). Note, that (1 + 6(r))/2 = (^)^. Thus we need to check that 

2( )^ - 1 < 1 - (— ^-r < = > 2(^— )?=3 < 2 - (— ^— ■ )?=3 

r — 1 r — 1 r — 1 r — 1 

<=> ( -V- 2 < 1 



r-1 ( r _ 1)1+^2 

< = > -)r-a < 

r — 1 r 

<=> 2< (r-l)(l-i) r - 2 . 

r 

Note that g(r) = (1 — ^) r ~ 2 is a strictly growing function of r since log(g(r)) = (r — 
2) log(l — i). Therefore, it suffices to verify that for r = 5 we have 2 < 4(1 — |) 3 « 2.048. 
□ 

Applying Lemma ^2 & = & is the unique solution of /(&) = 6. When r = 2, we obtain 
6=1/3. When r = 3 we obtain b = 1 - (1 + bf/A or b = 2VS - 3 » .4641 . When r = 4 
we find numerically that 6 « .5419. 
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Figure 2: /(/(&)) = b has more than one solution when r = 8. 



Recall, that the distribution F* of B is given by ¥(B = 0) = b and P(B > t) = (1 - 6)e~*. 
Applying of Theorem El and generating B\, . . . , B r ~ F* independently, we obtain 

(45) lim X w (n,r) = Z" 00 ^ > f 

We again condition on j out of r terms 5^ to be equal to zero and the rest positive. 
This occurs with probability ( r )^'(l - b) T -i . W.l.g. assume that the non-zero terms are 
£?x, . . . , B r _j. We have 

/co poo 
te- l ¥(t > B i) dt =1 ~ te ~* P ( Yl Bi > *) dt 

l<i<r-j l<«<r-i 

Repeating the calculations pS|). we obtain that P(Xa<?:<r-j ^ > ^) = Eo<i<r-i-j fr e ~*- 
Then the expression of the right-hand side of ()46|) becomes 



(47) I- I te-t V tr*d* = l- V l+i = r -^ + 2 . 

^ u 0<i<r-l-j 0<i<r-l-j 

Combining, we obtain 

iim x - (n ' r) = y ( r )v(i-br jr ~ j \ 2 = -2 y f r V'( 1 --r-+ 

0<j<r VJ/ 0<j<r VJ/ 



0<j<r VJ/ 



22 



The first summand is simply (\ + \) r ■ We compute the second summand using the following 
probabilistic argument. Rewrite the expression as 

0<j<r V/ ^2 T 2> 

The sum above is simply the number of successes in r Bernoulli trials with the probability 
of success equal to (± - §)/(± + §) = (1 - 6)/(l + 6). Namely, it is r(l - 6)/(l + 6). We 
obtain 

(48) lim J - (n ' r) = (I + + = (1 + &) r 1 (2 + 26 + r - r6). 

v ; n n v 2 2 J K 2 + 2b' 2 r + 1 v ; 

Recall from Lemma that (1 + b) r ~ 1 /2 r ~ 1 = 1 — 6. This proves (fT2*|) Plugging the 
corresponding value of 6 for r = 2, 3, 4, we obtain lim n Jw< ^' 2 ^ = 2 lim n Iw (™> 3 ^ = .6077 and 
lim n Jl " ( ' n ' 4 ^ = .5632. This concludes the proofs for the case of maximum weight independent 
set in Gzin) and G 4,(71). Before we continue with other cases, it is convenient to prove 
Corollary HJ as the proof is almost immediate from above. 

Proof : [Proof of Corollary 0] We have established above that Tj r has a unique fixed 
point iff r < 4. Applying ^TUJ) we need to compute E[1{W - Y.i<i< r B i > °}] = > 
J2i<i< r Bi), where W ~ Exp(l), Bi ~ F*. Instead of computing this quantity directly, note 
that the probability above is exactly 1 — 6, if we the summation above was up to r — 1 not 
r. Repeating the computations up to ()39|). we obtain 

f(w >j2 B ^ = 

l<i<r 

Plugging the obtained values of 6 for r = 2, 3, 4, we obtain (|18|) . □ 
Independent sets in G(n,c/n). 

Let F* denote now any fixed point of T 2 = T| c . We introduce again 6 = F*(0) = 

F(W - J2i<i< m Bi ^ °)> where ^ ~ T ( F *) and m ~ Pois(c). Similarly, 6 = T(F*)(0) = 
F(W — J2i<i< m Bi < 0), where 5j ~ F*. Repeating the computations done for G r (n), we 
obtain similarly to ()39jl that conditioning on m = k 

6 = 1- ( r. 

v 2 ; 

Then 6 = £ fc > ^e~ c (l - (±±^) fc ) = 1 - e"^. Similarly, 6=1- . Thus, 6 

must satisfy 1 — 6 = exp(— | exp(— 1(1 — 6))). Recall, however, that by the second part of 
Theorem ^ the equation above has a unique solution iff c < 2e. In this case 6 is also the 
unique solution of 1 — 6 = exp(— 1(1 — 6)). We now apply (J7J) of Theorem |2] to compute 
lim n Xw \ n ' c ) ; where we substitute m ~ Pois(c) for r. In order to shortcut the computations, 
we use flH3) and We have 

^5^0 r te - P(i> ^ Bi)rft! 

fc>0 ' 70 Ki<fc 
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where Bi ~ F*. Recall, though, from (|48j). that 

fc(l-fe) 



fOC 1 _|_ 7 

/ te-*p(f > X] = (^) fc (l + 



2(1 + b) 

Ki<k 



)• 



Combining 



,. X w (n,c) c(i-i>) (1 — 6) c(l + 6) c(i-i,) c(l — 6). <=(i-io 

hm v ; = e ^ + -+ -e ~ = (1 + — -)e ~. 

n 2(1 + 6) 2 v A ' 

Using 6=1— exp(— c( - 1 ~^ ), we obtain (JHJ). 
Matchings in G r (n). 

Let .F* denote any fixed point distribution of T = T^ if . Then we have the following dis- 
tributional identities B = maxi<;< r _i(0, Wi - Bi), Wi ~ Exp(l), B { ~ F*, 5 ~ T(F*) and 
B = max^^^O, - A), ~ Exp(l), B { ~ T(F*), B ~ F*. Let S = F(W - B < 0), 
where VP ~ Exp(l) and B ~ T(F*), and let 6 = P(W - B < 0), where 5 ~ F*. Then for 
any t > 

(49) P(£<t)=P( max (W* - < t) = (1 -P(W X > t + B 1 )) r - 1 = (1 - e _ *(l - 6)) r_1 . 

l<i<) — 1 

In particular 

(50) P(5 = 0) = 6 r " 1 , dF{B <t) = {r-l){l-b)e- t {l-e- t {l-b)) r - 2 dt, t>0 
(We note as above that F(W\ = ■) = 0, since W\ has a continuous distribution). Then 

poo 

b = e- l ¥(B > t)dt 
Jo 

e-*(l-(l-e-'(l-6)) r - 1 )dt 

= 1- / e -*(l - e _t (l - tyy^dt 
Jo 

poo 

= 1- / e-'(l -e-'(l -6)) r - 1 rft 
Jo 

POO 

= 1- / (l-e-'(l-6)) r -M(-e-*) 

= 1- [ {1- z{l-b)) r ~ l dz 
Jo 

l (i - z(i - b)y 
K-(i-S)) 

1 - 6 r 

(51) = ] 



r(l - 6) 
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Similarly, we obtain 



r(l-6)' 

and combining we conclude that b must be a fixed point of the equation f(f(x)) = x, where 

m = i - — - At. 

r(l — a; J 



Lemma 12 For ei>en/ r > 2 i/ie equation f(f(x)) = x has a unique solution x* in the 
range [0, 1], which is the unique solution of the equation f(x) = x. 

Proof : Note that (1 — x r ) /(l — x) — 1 + x + • • • + x r ~ 1 and therefore f(x) is a strictly 
decreasing function with /(0) = 1 — 1/r, /(0) = 0. Therefore /(#) = x has exactly one 
solution x*. We now prove that f(f(x)) > x for all x < x* and f(f(x)) < x for all x > x*. 
This would complete the proof of the lemma. We need to show that for x < x* 

1 _ fr( x ) 

/(/(*)) = 1 - J f J > X, 

which is equivalent to 1 — f r (x) < r(l—f(x))(l—x). But since / is a decreasing function and 
f(x*) = x*, then f(x) > x for all x < x* and therefore 1— f r (x) < l—x r = r(l—f(x))(l—x). 
Similarly, when x > x* we have f(x) < x* and then 1 — f r (x) > l — x r = r(l — f(x))(l— x), 
resulting in 

/(/(*)) = i - )~ f y\ < x. 

r (! - f{x)) 

□ 

We conclude that b = b is determined as the unique solution of 

1 - b r 

(52) b = 1 



r(l-b) 



and the unique fixed point of T 2 is the distribution given by (J49j) with b = b given above. 
Now, using (JHJ) of Theorem |21 and (JoTty we have 



lim M ^ r ) = E\Y Wil{Wi - Bi = max (0, W< - BA}}, 

n n Kj<r 

l<i<r 

= rE[W^l{Wi - B 1 > max (0, W 3 - BA}], 

2<j<r 



(53) = r 



POO 

/ te-*6 r_1 P(t > max(Wj - Bj))dt 

JO 2<j<r 



pco rt 

+ r / te~\r- 1)(1 — 6) x 
(54) x e~ z (l - e~ 2 (l - b)) r ~ 2 F(t -z> max (^ - Bj))dtdz, 

2<j<r 
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where the summands and (J53j) corresponds to conditioning on B\ = and B\ = z > 0. 
We now compute the integrals in these summands. We have 

poo poo 

/ te- t 6 r_1 P(t > max(Wj - Bj))dt = / te^V'^- 1 (t > W 2 - B 2 )dt 
Jo 2 ^ r Jo 

POO 

= / te-*6 r - 1 (l - P(W 2 > t + B 2 )) r - l dt 
Jo 

POO 

= / te-% r -\l - e- t ¥(W 2 > B 2 )) r ~ l dt 
Jo 



OO 



te-*6 r - 1 (l-e-*(l-6))) r - 1 dt 
Similarly, we obtain that the integral in (JHlj) is equal to 

*e"*(r - 1)(1 - b)e- z (l - e~ z (l - &)) r ~ 2 (l - e _t (l - b^Y^dtdz 
te~\r - 1)(1 - 6)(1 - e- z (l - &)) r ~ 2 (l - e _t (l - b))) r - l dtd{-e- z ) 




'0 ./o 

"CO 




./0 



te _t (r - 1)(1 - 6)(1 - e _t (l - 6)))'- 1 cft / (1 - iw(l - &)) r ~ 2 du> 

I! JO 

te-*(l - e~*(l - fo))) r - 1 dt((l - iu(l - ft))'- 1 ° ) 



o 



oo />oo 



te"'(l - e-*(l - ft))) 7 "- 1 ^ - / te-*(l - e -*(l - b))) 2r ~ 2 dt 
o Jo 

Substituting to summands in (|53|) and (|51|i. we obtain (JT5j) . 
Matchings in G(n,c/n). 

The derivation is very similar to the one for A4 w (n, r). We introduce 6 and 6 exactly as 
above. The equation (JoT]) becomes 



^o m! (m + l)(l-6) ; 



c »n+l j ^ c »n+l^m+l 



1 n m+i i 

i V — e" c + - V 



c(l-&)fTo^ m+1 ) ! c(l-S)^ {m+l)\ 

1 e -c(l-6) 
= 1 ^ + -, 

c(l - b) c(l - 6) 

Then 1 — b is a fixed point of the equation f(f(x)) = x, where f(x) = x ~ e cx cx ■ That is 

l-e~ cx 

i ^ c c 

1 —p-cx 

c— - — 
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First we note that x = does not satisfy f(f(x)) = x, so we assume x > 0. The expression 
above then becomes e~ cx + ex 2 — 1 = 0. The function e~ cx + cx 2 — 1 is a strictly convex 
function which is equal to zero at x — 0, has derivative — c < at x = and diverges to 
infinity as x diverges to infinity Therefore its graph has exactly one intersection with the 
horizontal axis in x G (0, oo). Note that at x — 1 the value is e~ c + c — 1 > (for every 
c > 0) therefore there exists exactly solution to e~ cx + cx 2 — 1 = in x G (0, 1). We conclude 
that b is uniquely determined by the equation e - ^ 1- ^ + cb 2 — 1 = 0. The remainder of 
the calculation is done just like for G r (n) by conditioning first on specific values of r and 
recalling m = r with probability y!e~ c . We omit the fairly straightforward calculations. □ 

5.3 Extensions: Cycles and (non-random) low degree graphs 

In this section we first prove CorollaryEl We then focus on proving TheoremHJ thereby improving 
the previous bound of Hopkins and Staton [HS82J. 

Proof : [Proof of Corollary E] The result is essentially established en route of proving 
Theorems 121 and El Recall that in the proof of Theorem El the only place we used the randomness 
of our regular graph G r (n) was to say that a constant depth-ci neighborhood H(d) of a randomly 
selected node i G [n] is a depth-ci r-regular tree w.h.p. In the case of a cycle any constant depth-ci 
neighborhood of any node is w.p.l a path of length 2ci with the selected node in the middle. This 
is depth-ci 2-regular tree. The answer for the maximum weight matching on a cycle is the same 
as for independent set since we may simply assume the weights are assigned to nodes which are 
to the left for each edge. □ 

Proof : [Proof of Theorem 0] As in the proof of Corollary El we recall that the only place 
the randomness of G r (n) was used in the proof of Theorem [2] and Theorem El was the fact that a 
constant depth-ci neighborhood of a randomly chosen node i G [n] is w.h.p. a depth d r-regular 
tree. In order to continue the proof we introduce a class of almost r-regular (non-random) graphs 
with large girth. Given positive integers d, R let Q r {n, d, R) denote the class of n-node graphs with 
girth at least d and such that all but at most R nodes have degree r and the remaining nodes have 
degree less than r. Given G G Q r {n, d, R) a random node % G [n] its depth-ci neighborhood H(i, d) 
is w.h.p. (with respect to randomness of the choice of i) a depth-ci r-regular tree. Repeating the 
proofs of Theorem |21 and El and Corollary 0] we conclude that for the case r = 3 for every constant 
R and for all sufficiently large n, ci, the expected maximum weight of an independent set in G is 
given by (fT3*|) with r = 3, that is ~ .6077n. Also by Corollary 01 the expected cardinality of the 
independent set achieving this weight is given by (|18|) with r = 3, that is ~ .3923n, also for all 
n and ci sufficiently large. Then the maximum cardinality of an independent set in G is at least 
.3923rt — o(n) and this is a non-probabilistic bound since we consider a maximum cardinality 
independent set in a fixed graph G. Thus for every constant R 

KG) 

liminf liminf min > .3923.. . 

d n Geg 3 (n,d,R) 71 

To finish the proof we need to obtain similar lower bound for the class Q(n, 3,ci) - all graphs 
with degree at most r = 3 and girth at least ci. Suppose we are given any such n-node graph 
G. Consider any two nodes which have degree < 3 such that i 2 ^ H(ii,d) (if any such 
pair exist). Connect by an edge. This operation can only decrease the size of the largest 
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independent set in G. We claim that in addition the resulting graph still has girth at least d. 
In fact, if the added edge becomes a part of a cycle with length less than d, then in the original 
graph the nodes i\, i 2 were connected by path of length at most d — 2, contradicting the condition 
i 2 ^ H(ii, d). Therefore the resulting graph still has girth at least d. We continue this operation 
for every such pair i\, i 2 until we either do not have any nodes with degree less than 3 at all or 
all such nodes belong to H(i, d) for some node i. Since \H(i, d)\ < 1 + 3 + 3 2 + . . . 3 d = R, we 
obtain a graph G' € Q r (n, d, R). Then 

1{G) 1(G) 
liminf liminf min > liminf liminf min > .3923n. 

d n G&g(n,3,d) 71 d n Ge5 3 (n,d,K) U 

This completes the proof of the theorem. □ 



5.4 Deterministic and Bernoulli weights 

Are the results obtained above relevant to the case when the weight of each node and edge W — 1 
with probability one or, in general, when the weights take some discrete values? Let us examine 
these questions with respect to our usual four operators T. We start with the case W — 1. For 
T = Tj r , the corresponding distributional equation is B = max(0, 1 — J2i<i< r -i ^»)- ^ ^ = 0> 
w.p.l, then B — 1, w.p.l, and vice verse. Thus B = and B = 1 are two fixed points of T 2 , and 
T 2 does not have a unique fixed point distribution. There is a physical explanation for the lack of 
uniqueness, coming directly from the lack of long-range independence. Given a depth-c? r-regular 
tree T with all the weights equal to unity, note that the boundary does carry a non-vanishing 
information about the root in the following sense. If all the leaves of the tree (boundary nodes) 
are conditioned to belong to the maximum weight independent set, then all the parents of leaves 
cannot be part of the set. Then the maximum independent set is obtained by selecting all the 
nodes in level d, not selecting nodes in level d — 1, selecting all the nodes in level d — 2, and so on. 
In the end whether the root is selected is fully determined by the parity of d. Thereby we do not 
have a long-range independence. Contrast this with the discussion in Brightwell and Winkler 
BW03J, where similar observation is used to show long-range dependence for Gibbs measures 
on infinite regular trees for the hard-core model. It is not hard to see that the similar lack of 
long-range independence holds for maximum weight matchings, where the weights are all 1. 

The situation is different for T = Tj jC . Let F be a distribution function given by F(t) = 
p, t e [0, 1),F(1) = 1. Namely, F is simply a Bernoulli distribution with parameter p (Be(p)). 
If B = max(0, 1 - J2i<i< m B i)> where B { ~ F and m ~ Pois(c), then B = 1 if J2i<i< m B i = °> 
which occurs with probability ^ k>0 jr e ~ c P k — e~ c<yl ~ p \ and B = otherwise. Thus T(F) is 
Be(pi) where pi = 1 — e~ c ^~ p \ Similarly T 2 {F) is Be(p2)> with with p 2 = 1 — e~ ce cl P) . In 
general, for s = 1,2,..., T 2s (F) is Be(p 2 s) with 1 — p 2s = e~ ce c<1 P2a ~ 2 \ By Proposition [T] we 
know that for F = Be(0) and F = Be(l), T 2s (F) converges to some fixed point distributions 
F**,F** which by the argument above are Be(p**), Be(p**) with both p = p^ and p = p** 
satisfying 1 — p = e~ ce P . Recall from Theorem ^ that the equation x = e~ ce cx has a unique 
solution iff c < e. By Proposition ^ this implies that when c < e, T s (Fq) converges to Be(p*) for 
any starting distribution Fq. Again applying Proposition [T] we obtain that T S (F ) converges to 
Be(p*), T and T 2 have the same unique fixed point distribution - Be(p*), where p* is also the 
unique solution of 1 — p* = e~ c ^ 1 ~ p *\ 
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It is a simple exercise to see that the same holds for T = Tm,c' T 2 has a unique fixed point iff 
c < e, in which case the fixed point distribution is also Be(p*). This is, of course, fully consistent 
with Theorem ^ We summarize these observations. 

Proposition 2 Let the nodes and edges weights be equal to one, w.p.l. For every r > 1, 
T| r , Tjb r have at least two fixed point distributions. T| c , c have the unique fixed point distri- 
bution iff c < e in which case the unique fixed point distribution is Be(p*), where p* is the unique 
solution of 1 — p* = e~ c ^~ p *\ 

Can we fully reproduce Theorem^ for the case c < e? The problem with the case F w — 1 as, 
generally, with non-continuous distributions F w , is that the probability of W— J2i<i< r = is no 
longer zero. As a result we do not have an exact condition purely in terms of B = max(0, W — 
J2i<i< r Bi) for whether the root node belongs to say maximum weight independent set. One 
natural approach would be to approximate F w with a continuous distribution. But the difficulty 
is the lack of closed form expression for the solution of fixed point of T 2 (F*) = F*. Such a solution 
F* can, though, also be approximated by computing T 2s (0) and T 2s (F w ) (T 2s (F w ^ r ),T 2s (F w , c )) 
for matching) for large s such that the differences between the two distributions is sufficiently 
small. 

Suppose now the weights W{ of the nodes are distributed as Be(z) for some parameter z G 
[0, 1]. For simplicity we will only consider the case r = 3 and T = Tj 3 and obtain a complete 
criteria for uniqueness. 

Proposition 3 For r = 3 the operator T| r _ 1 has a unique fixed point distribution iff z G [|, 1]. 
In this case the fixed point distribution is Be(p) with p = ^("l-z") 1 1 

Proof : Let F = Be(p) for any x G [0, 1]. Then for B' ~ T(F) we have B' = max(0, W - 
B\ — B-i) where -Bi,i?2 are independent and distributed as F. Then B' = when W = or 
when W — 1 and B\ + B% > 0. This occurs with probability z + (1 — z)(l — p 2 ) = 1 — (1 — z)p 2 . 
Thus B' ~ Be(l — (1 — z)p 2 ). Repeating the development above for the case of the deterministic 
weight, we need to analyze the number of the solutions to the equation f(f(x)) = x where 
f(x) = 1 — (1 — z)x 2 . First, the equation f(x) = x leads to the unique solution x* = ~^0^r^ 
(it is simple to check that the solution is in [0, 1] for every z G [0, 1]). Also x* is a solution of 
g{x) = f(f(x)) = x. We need to show that this is the unique solution of this equation iff z > 1/4. 

First we show that when z < 1/4, g(x) > 1. Since g(0) > this would imply that there exists 
a solution of g(x) = x in the open interval (0, x*) and the case z < 1/4 would be resolved. We 
have g(x*) = 4(1 — z) 2 f(x*)x* = (2(1 — z)x*) 2 , where we use f(x*) = x*. Then the inequality 
g(x) > 1 holds when x* > 1/(2 — 2z) which after as simple algebra is reduced to z < 1/4. Thus 
when z < 1/4 there are more than one fixed points of T 2 = T| 2 . 

Suppose now z > 1/4. Let again x* be the unique solution to f(x*) = x* and consider 
g(x) = 4(1 — z) 2 f(x)x. We claim that when z > 1/4, g(x) < 1 for all x G [0,1] and when 
z = 1/4, g(x) < 1 for all x ^ x* and g(x*) = 1. This immediately implies that g(x) = x has at 
most one solution, meaning it has exactly one since g(x*) = x*. Note g(0) = < 1. We now 
prove that g(l) = 4(1 — z) 2 z < 1 and g(x) < 1 for all x such that d ^ — 0, except for x = x* 
for which we will show that g(x*) = 1. 
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We start with g(l) = 4(1 — z) 2 z = <p{z). Note 0(0) = 0(1) = < 1. Then the maximum 
is achieved at the points z where <j>(z) = —8(1 — z)z + 4(1 — z) 2 = 4(1 — z)(l — 3z) = 0. We 
have already considered the case z — 1. Otherwise z = 1/3 for which 0(1/3) = 16/27 < 1. Thus 
g(l) < sup 2e[0>1] (j){z) < 1. 

Consider now points x such that = 4(1 — z) 2 {— 2(1 — z)£ 2 + 1 — (1 — 2;)x 2 ] = 0, from 

which we obtain x = l/y/3 — 3z. Plugging this into 



g(x) = 4(1 - z) 2 (l - (1 - z)x 2 )x = 4(1 - z) 2 - 



2 1 64(1 -z) 3 , 

3 73^3^ ~ ^ 27 ' 



The last expression is smaller than unity whenever z > 1/4. When z = 1/4 the expression is 
equal to the unity. In this case, however, x = 1/V3 — 3z = 2/3 which is checked to be equal to 
x* = ^2(1- z) 4 * w hen z = 1/4. This concludes the proof of the proposition. □ 

6 Conclusions 

We have derived in this paper the limits of maximum weight independent sets and matchings 
in sparse random graphs for some types of i.i.d weight distributions. Our method is based on a 
certain local optimality property which states loosely that in cases of certain distributions of the 
random weights, the optimal random combinatorial structure under the consideration exhibits 
a long-range independence and, as a result, the value which each node (edge) "contributes" to 
the optimal structure is almost completely determined by the constant depth neighborhood of 
the node (edge). We certainly believe that such local optimality holds for many other random 
combinatorial structures and it seems to be an interesting property to study by itself, not to 
mention its applications to studying random combinatorial structures. 
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